Exact conditions in finite temperature density functional theory 
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Density functional theory (DFT) for electrons at finite temperature is increasingly important in 
condensed matter and chemistry. Exact conditions that have proven crucial in constraining and 
constructing accurate approximations for ground-state DFT are generalized to finite temperature, 
including the adiabatic connection formula. We discuss consequences for functional construction. 



Because of the small mass ratio between electrons and 
nuclei, standard electronic structure calculations treat 
the former as being in their ground state, but routinely 
account for the finite temperature of the latter, as in ab 
initio molecular dynamics [1]. But as electronic struc- 
ture methods arc applied in ever more esoteric areas, the 
need to account for the finite temperature of electrons in- 
creases. Phenomena where such effects play a role include 
rapid heating of solids via strong laser fields [2] , dynamo 
effects in giant planets [3], magnetic [4, 5] and supercon- 
ducting phase transitions [6, 7], shock waves [8, 9], warm 
dense matter [10], and hot plasmas [11-13]. 

Within density functional theory, the natural frame- 
work for treating such effects was created by Mermin [14] . 
Application of that work to the Kohn-Sham (KS) scheme 
at finite temperature also yields a natural approximation: 
treat KS electrons at finite temperature but use ground- 
state exchange-correlation (XC) functionals. This works 
well in recent calculations [8, 10], where inclusion of such 
effects is crucial for accurate prediction. This assumes 
that finite-temperature effects on exchange-correlation 
are negligible relative to the KS contributions, which may 
not always be true. 

The uniform electron gas at finite temperature (also 
called the one-component plasma) has been well-studied, 
and has in the past provided the natural starting point 
for DFT studies of such finite-temperature XC effects, 
as input into the local density approximation (LDA) at 
finite T [15]. However, the LDA is too inaccurate for 
most modern applications of DFT, and almost all recent 
calculations use a generalized gradient approximation or 
hybrid with exchange [16]. The errors of LDA would 
typically be enormous relative to the temperature correc- 
tions we seek, especially for correlation, and so could lead 
to quite misleading results. Accurate calculation of finite 
temperature contributions requires accurate approximate 
functionals. Magnetic phase transitions bear an addi- 
tional difficulty: The low-lying excitations are collective, 
i.e., magnons whose description requires non-collinear 
version of spin-DFT. Hence, a finite-temperature version 
of spin-DFT involving only spin-up and spin-down den- 
sities and thus only spin-fiip excitations, is bound to fail 



in predicting, e.g., the critical temperature [4]. 

The most fundamental steps toward both understand- 
ing a functional and creating accurate approximations 
are deriving its inequalities from the variational defini- 
tion of the functional. These yield both the signs of en- 
ergy contributions and, via uniform scaling of the spatial 
coordinates, basic equalities and inequalities that non- 
empirical functionals should satisfy by construction. The 
adiabatic connection formula [17] is intimately related. 
Here, we (i) establish the fundamental functionals needed 
for treating finite temperature, (ii) prove the most basic 
properties (signs of the energy contributions), (iii) show 
that the temperature must be scaled simultaneously with 
the spatial coordinate, (iv) derive the inequalities under 
such scaling, and (v) give the adiabatic connection for- 
mula for finite temperature. These results establish the 
basic rules for all finite-temperature KS treatments. 

Central to the thermodynamic description of many- 
electron systems is the grand-canonical potential, defined 
as the statistical average of the grand-canonical operator 



n = H - tS - ixN, 



(1) 



where H , S, N, t and ^ are the Hamiltonian, entropy, 
and particle- number operators, temperature and chemi- 
cal potential, respectively. In detail, H = T -\- T4e + V, 
where T and Vee are the kinetic energy and the Coulomb 
electron-electron interaction operators, and V represents 
an external scalar potential v{t). The entropy operator 
is given by S* = — A: In F , where k is the Boltzmann con- 
stant and F = ^ wnaY^ N,i) (^iv.il is a statistical oper- 
ator, with \^ N,i) and WN,i being orthonormal A^-particle 
states and statistical weights, respectively, with the latter 
satisfying the (normalization) condition '^jq .iWN,i = 1- 
The statistical average of an operator A is obtained as 



A[t] = Tr {f i} = Y,Y. wnA'^nAM'^n. 



(2) 



N 



The thermodynamical equilibrium properties of many- 
electron systems are obtained from the knowledge 
of the grand-canonical statistical operator F^ = 



E 



N,i '^N.i 







in 



)(^'^ J, where |^'^^.j) are the TV-particle 
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eigenstates of H with energies ^ , and the equi- 
librium statistical weights are given by ,- = 

^'^'^P'"^^^-,:^^^' where /? = f [18]. The Gibbs 

principle ensures that r° minimizes the statistical aver- 
age of the grand-potential operator. We emphasize that 
r° is unique [14] and that in the limit of zero tempera- 
ture, for systems with degenerate ground states, it leads 
to ensembles with equal statistical weights. 

To create a DFT at finite temperature, Mermin [14] 
rewrites this as (in modern parlance) 



a;_^=mm|F-[7i] + j 



(Pr n{r) {v{r) — yu) 



(3) 



where the minimizing n{r) is the equilibrium density 
n°(r), and 

F^[n] := min i^^[f ] = min |r[f] + Vcc[t] - TS[t]] , 



(4) 

is the finite-temperature analog of the universal 
Hohenberg-Kohn functional, defined through a con- 
strained search [18, 19]. This depends only on r and 
not on fi. We denote r'^[n] as the minimizing statistical 
operator in Eq. (4), and define the density functionals: 

S-[n] S[r^[n]], (5) 

i.e., each density functional is the trace of its operator 
over the minimizing T for the given r and density. 

Next consider a system of non-interacting electrons at 
the same temperature r, and denote its one-body poten- 
tial as Vs{r). All the previous arguments apply, and we 
choose Vs{r) to make its density match that of the in- 
teracting problem. This defines the KS system at finite 
temperature. Because it arises so often in this work, we 
define the kentropy as 



K^[r] :=r[f] -TS'[f] 



(6) 



and we show it plays an analogous role to the kinetic 
energy in ground-state DFT, to which it reduces as t — 
0. The non-interacting functional is just 



(7) 



from Eq. (4) applied without Vcc, and we define: 

T,nn]:^T[ri[n]], S;[n] := S[ri[n]]. (8) 

Next we define the difference functionals that are crucial 
to the KS method. Write 



(9) 



where Lf^ [n] in terms of the density has the form of the 
usual Hartree energy, and expressing rixl"-] terms of 



the module square of the one-body density matrix stem- 
ming from ri[n,] [20] we observe that ^'^[n] < 0. 
The kinetic correlation is 



T^[n] ■.^T[t-[n]]-T[tl[n] 



(10) 



and similarly define [n] and [n] , while the potential 
contribution is 



U^[n] :=Kc[r"[n]]-1/ec[r;[n]: 



(11) 



The sum of the energy components is, as in ground-state 
DFT, the correlation energy, E^[n] T^[n] + U^[7i], 
while the grand-canonical correlation potential is 



(12) 



and 5^xc["' 



We now prove the most basic theorems about the signs 
of our quantities. To show that the correlation-kentropy 
(or kentropic correlation) is always positive, begin by 
noting i^'^[r3[n]] < if'^[r'^[n]], because rs[n] minimizes 
i4r'^[r]. By inserting the definition, Eq. (6), we find 
^^['T'] ^ 0, with equality only when the interaction is 
zero. It is the kentropic correlation that is guaranteed 
to be positive, not the kinetic correlation alone, con- 
trary to the pure ground-state case [21]. Similarly, since 
F^[t^[n]] < F^[f^[n]], we find ni[n] < 0. Combining 
these results with Eq. (12) implies U^[n] < 0. Thus 

ni[n] < 0, < 0, U^in] < 0, K^in] > 0, (13) 

and no approximation should violate these basic rules. 

Some of the most important results in ground-state 
DFT come from imiform scaling of the coordinates [21, 
22]. In the following considerations, when we refer ex- 
plicitly to wavefunctions, we shall restrict to wavefunc- 
tions having finite norm on their entire domain of def- 
inition. Under norm-preserving homogeneous scaling of 
the coordinate r — 7r, with 7 > 0, to the scaled wave 
function [21] 



1'''(ri,...,rjv) := 7^'^*(7ri, 7i-Ar), 



(14) 



corresponds the scaled density ny{r) = ^^n{'~fr). Writ- 
ing \I''''(ri, Tat) = (ri, rAr|\E'''') in terms of the 
(representation-free) element 1^'') of Hilbert space, the 
scaled statistical operator is defined as 



EE 



WN. 



N,: 



(15) 



N 



where the statistical weights are hold fixed, i.e., the scal- 
ing only acts on the states. 

With the above definition, the statistical average of an 
operator whose pure-state expectation value scales ho- 
mogeneously [21], scales homogeneously as well. In par- 
ticular, we have: T[f^] = 7'r[f], V,,[t^] = 7Kc[r], 
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N[t^] = iV[f], and S[t^] = S[r]. The scaling behavior 
of the density functionals is, however, more subtle. First 
consider the non-interacting functionals in some detail. 
Because Tl[n] minimizes K'^ , Eq. (7), and 

A-[f,] = 7' (rit] - S[t]^ = 7'^^/"'[f] , (16) 

then 

tl[n,] = tl/f[n] , F;K] = 7' ^^s^/"' W • (17) 

In particular we notice that 

S:[n,]=S:^'^'[n]. (18) 

For non- interacting electrons, the statistical operator at 
a given temperature that is the minimizer for a given 
scaled density is simply the scaled statistical operator, 
but at a scaled temperature, an effect that is obviously 
absent in the ground-state theory. 

There arc further simple implications. First, if we in- 
vert the sense of Eq. (17), we can write: 

F^'[n]^^F,nn^], (19) 

i.e., knowledge of at any one finite r generates its 

entire temperature dependence, via scaling. Furthermore, 
it must always collapse to the ground-state KS kinetic 
energy under scaling to the high-density limit: 

Ts[n] = lim F;K]/72. (20) 

Similarly, in the low-density limit 

5rW = -limF;K]/r, (21) 

7—^0 

where S'^[?t.] is the non- interacting KS entropy in the 
high-temperature limit. 

Next, we consider the interacting case. The exchange 
contribution is much simpler than correlation, because 
it is extracted from the one-particle density matrix. Be- 
cause F.^ [r] and U[r] scale linearly with 7, and using the 
simple scaling relation for Fg, Eq. (17), 

ni[n^]=jnl/^\n]. (22) 

This scaling result is important in ground-state DFT, 
where it restricts the dependence of the exchange- 
enhancement factor to depending on just the reduced 
density gradient [22]. But the more interesting case is 
correlation. From the definition, Eq. (4), 

F-[n,]<F-[r;'[n]], (23) 

since Tl^ [n] has density n-y, and r' is any temperature. 
Using the scaling properties and choosing r' = T/7^, then 
the fundamental inequality of scaling is 

A- [n^] + K] < 7' i^^/-^' N + 7 Ko^'^' N • (24) 



To find a condition on the kentropy alone, define n'(r) = 
n^(r), 7' = 1/7, and r' = t/j'^ in Eq. (24). Multiply 
the result by 7', and combine with Eq. (24), to find 

K'' K] < 7^ if^/'^' [n] , 7 > 1 . (25) 

This is the finite temperature analog of the subquadratic 
scaling of the kinetic energy in the real system [21]. An- 
other combination isolates the repulsive contributions: 

VJM'>]>lVjJ'^"[n], 7>1. (26) 

These inequalities loosely constrain the behavior of these 
large energies. Much more important is to subtract out 
KS quantities that scale simply, to find for 7 > 1: 

K^c K] < 7' K'J''^ [n] , U^a K] > 7 U^^''^ [n] . (27) 
One more application of Eq. (24) yields 

ni[n^]>jnl/''\n], 7>1, (28) 

the fundamental scaling inequality for the correlation 
contribution to the grand canonical potential. The in- 
equalities, Eqs. (25-28), which are reversed if 7 < 1, pro- 
vide tight constraints on these functionals and arc rou- 
tinely used in non-empirical functional construction in 
the ground state[22]. For example, combining Eq. (22) 
with Eq. (28) in the high-density limit, yields: 

ni[n] = lim nZlr[ny]/j. (29) 

7—^00 

This scaling procedure can usually be applied easily to 
any approximate fi^c [n] to extract its separate exchange 
and correlation contributions. 

Lastly, we consider the adiabatic coupling constant for 
finite temperature, its relationship to scaling, and derive 
the adiabatic connection formula. Define 

F^'^n] = min |T[f] + A Kc[f] - TS[t]} , (30) 

with r'^'^[n] being the corresponding minimizing F. By 
scaling, it is straightforward to show: 

f--^ [n] - f [m/A] , F^'^[n] = F^/^' [^/a] • 

(31) 

where quantities with one superscript are evaluated at 
A = 1. Eq. (31) is the interacting generalization of Eq. 
(17) and shows that, even in the presence of interactions, 
simple equalities are possible, but at the price of altering 
the coupling constant. In particular we notice that 

S^^^[n] = S^/^'[n,^^]. (32) 

Of course, non-interacting functionals are not affected by 
a coupling constant modification. Eq. (22) implies that 
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the exchange and Hartree density functionals have a hn- 
ear dependence on A. Employing the minimization prop- 
erty of Eq. (30) and the Hellmann-Feynman theorem, we 
find 



where 



(33) 



(34) 



Eq. (33) is the finite-temperature adiabatic connection 
formula, whose zero-temperature limit played a cen- 
tral role in ground-state DPT. U^^[n]{0) = ni[n] < 
(Eq.(13)), and the scaling inequalities can be combined, 
analogously to Ref. [21], to show that U^^,[n]{X) is mono- 
tonically decreasing in A. 

So far, all results presented have been exact. To see 
them in practice, consider the finite-temperature local 
density approximation (LDA) to r2Jc[?i] 



LDArr 



(35) 



ifr 



(n) is the XC grand canonical potential den- 
sity of a uniform electron gas of density n. Because a 
uniform electron gas is a quantum mechanical system, 
its energies satisfy all our conditions, guaranteeing by 
construction that LDA satisfies all the exact conditions 
listed here. In the Jacob's ladder of functional construc- 
tion [16], more sophisticated approximations should also 
satisfy these conditions. To give one simple example, Eq. 
(22) implies 



Lifr 



(n(r)) 



f(n(r))Fx(f(r)) 



where e™'f(n(r)) = -A^ n'^/^(r), = (3/47r)(37r2)i/3^ 
and f(r) = t jr^l'^iy:') is a dimensionless measure of the 
local temperature. Thus the largest fractional deviations 
from ground-state results should occur (in LDA) in re- 
gions of lowest density, but these contribute less in ab- 
solute terms. For a generalized gradient approximation 
(GGA), Eq. (22) implies 

a;««^-(n(r),|Vn|(r)) =er'(n(r))J^x(.(r),f(r)), (37) 

where the dimensionless gradient s is |V?7 |/(2fci?n) and 
kp = (37r^n)^/'^, i.e., the exchange enhancement factor 
Fx{s, f) depends on the temperature only via f. 

In summary, there is a present lack of approximate 
density functionals for finite temperature. We have de- 
rived many basic relations needed to construct such ap- 
proximations, and expect future approximations to either 
build these in, or be tested against them. In principle, 
such approximations should already be implemented in 
high-temperature DFT calculations, at least at the LDA 
level, as a check that XC corrections due to finite tem- 
perature do not alter calculated results. If they do, then 



more accurate approximations than LDA will be needed 
to account for them. 
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